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We present next-to-leading order (NLO) perturbative-QCD calculations of the cross sections for IN —> hX 
and IN —> jet Jf. The main feature of these processes is that the scattered lepton is not observed, so that the hard 
scale that makes them perturbative is set by the transverse momentum of the hadron or jet. Kinematically, the 
two processes thus become direct analogs of single-inclusive production in hadronic collisions which, as has 
been pointed out in the literature, makes them promising tools for exploring transverse spin phenomena in QCD 
when the incident nucleon is transversely polarized. We find that the NLO corrections are sizable for the spin- 
averaged cross section. We also investigate in how far the scattering is dominated by the exchange of almost 
real (Weizsacker-Williams) photons. We present numerical estimates of the cross sections for present-day fixed 
target experiments and for a possible future electron ion collider. 
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I. INTRODUCTION 

There has been growing interest recently, both experimen¬ 
tally [1-4] and theoretically [5-10], in the processes fN —> hX 
and fN jet V, the single inclusive production of a hadron or 
jet at large transverse momentum in lepton-nucleon scattering. 
In contrast to the far more customary process fN f'hX [11], 
for fN —» hX the scattered lepton in the final state is not ob¬ 
served, so that the process is truly one-hadron (or one-jet) in¬ 
clusive. The reason for the interest in fN —» hX comes from 
the study of single transverse-spin phenomena in hadronic 
scattering processes. It is well known that large single-spin 
asymmetries have been observed [12] for the process pp^ ^ 
hX, where p* denotes a transversely polarized proton. To ex¬ 
plain the large size of the asymmetries, and their persistence 
all the way from fixed-target to collider energies, has posed a 
major challenge to theory. Although a lot has been learned, 
it is fair to say that a fully satisfactory understanding has yet 
to be obtained. Measurements of corresponding asymmetries 
in the kinematically equivalent, but much simpler, processes 
fN'^ —> hX, fN'^ jetV have the promise to shed new light 
on the mechanisms for single-spin asymmetries in QCD. First 
fairly precise experimental data for fN'^ hX have recently 
been released by the HERMES [2, 3] and Jefferson Lab Hall 
A [4] collaborations. 

We note that at first sight one might consider the related 
process —> f'X (which is just the standard inclusive deep- 

inelastic (DIS) process) to be equally suited for transverse- 
spin studies in lepton scattering. However, the analysis of the 
corresponding single-spin asymmetry is considerably more 
complex because higher order QED effects are required for 
the asymmetry to be non-vanishing [13-17]. In the same spirit 
as tN'^ hX, also the processes fN'^ hX [18] with longi¬ 
tudinal polarization of the lepton and fN A*V [19] with a 
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transversely polarized A hyperon have been considered in the 
literature recently. 

The proven method for analyzing single-inclusive pro¬ 
cesses such as pp —» hX or fN hX at large transverse 
momentum rests on QCD perturbation theory and collinear 
factorization. Eor single-transverse-spin observables, this in¬ 
volves a twist-3 formalism in terms of three-parton correlation 
functions of the nucleon or the fragmentation process [20-29]. 
Interestingly, the recent study [29] suggests that the twist-3 
fragmentation effects could be the dominant source of the ob¬ 
served large transverse-spin asymmetries in pp^ —> hX. An 
alternative approach for describing the single-spin asymmetry 
in inclusive hadron production in pp* —» hX was devised in 
the context of a “generalized” parton model in which the de¬ 
pendence of parton distributions and fragmentation functions 
on transverse momentum is kept [30-32]. Although no such 
factorization in transverse momentum is known to be valid 
for a single-inclusive cross section, the approach has enjoyed 
considerable phenomenological success. 

Both the collinear twist-3 approach and the generalized par- 
ton model have been used to obtain predictions for the spin 
asymmetry in fA* —> hX. In Ref. [7] a leading order (LO) 
twist-3 analysis has been presented in terms of parton cor¬ 
relation functions that were previously extracted from data 
for pp* —> hX. The results obtained in this way fail to de¬ 
scribe the HERMES data [2, 3] for the spin asymmetries in 
fN^ hX. A comparison of perturbative calculations to the 
corresponding JLab data [4] is not possible as the data are for 
hadrons with transverse momenta below 1 GeV. The LO gen¬ 
eralized parton model approach, on the other hand, appears to 
give results quite consistent with the HERMES data [8-10]. 

In our view it is premature to draw any conclusions from 
these findings at LO. Given the kinematics (and the preci¬ 
sion) of the present data, one may expect higher-order QCD 
corrections to the cross sections and the asymmetry to be 
important [7] for a meaningful comparison of data and the¬ 
ory. At least next-to-leading order (NLO) corrections should 
be included. We stress that the twist-3 formalism, although 
so far only developed to LO, offers a well-defined frame¬ 
work for a perturbative study of the transverse-spin asym- 
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metry in (N'^' hX. This is in contrast to the generalized 
parton model, for which there is likely no systematic way 
of going to higher orders in perturbation theory. That said, 
NLO calculations within the twist-3 formalism are technically 
very challenging, and only a few NLO calculations have been 
performed for the simpler Drell-Yan [33] and semi-inclusive 
deep-inelastic scattering (DIS) cases [34]. 

In the present paper, we take a first step toward an NLO 
calculation of the transverse-spin asymmetry for iN'^ —> hX 
by computing the NLO corrections to the spin-averaged cross 
section for the process, which constitutes the denominator of 
the spin asymmetry. We present analytical results for the NLO 
partonic cross sections. To our knowledge, despite the vast 
amount of work performed for lepton proton scattering in the 
literature (see, for example [35-41]), this calculation has not 
been presented so far. We also present similar NLO calcula¬ 
tions for the process iN jetY. We note that the process 
(N ^ jet Y has also been extensively studied in terms of the 
concept of “l-jettiness” [42, 43]. Here one additionally writes 
the cross section differential in a variable t\ that characterizes 
the hadronic final state that is not associated with the produced 
jet or the nucleon beam remnant. In Ref. [43] the full NLO 
corrections for the l-jettiness were computed, where a fully 
numerical approach was adopted. In principle, it should be 
possible to recover our NLO results by performing a (numer¬ 
ical) integration over t\ of the results of [43]. 

Because of the propagator of the exchanged photon, the 
cross section for (N^ hX will contain contributions for 
which the photon is almost on-shell. This is not yet the case 
at LO where the high transverse momentum of the produced 
hadron requires the photon to be highly virtual. Starting from 
NLO, however, it may happen that the incoming lepton ra¬ 
diates the photon almost collinearly. This may then be fol¬ 
lowed by a 2 ^ 2 scattering process of the photon with a 
parton in the nucleon, which is perfectly capable of produc¬ 
ing the hadron at high P/,^. In processes where the scattered 
lepton is observed, such as iN i'hX, one can in fact se¬ 
lect such contributions by requiring the scattered lepton to 
have a low scattering angle. The incoming lepton then effec¬ 
tively acts merely as a source of quasi-real photons, and the 
process may be very accurately described in terms of a (per¬ 
turbative) distribution function for photons in leptons known 
as the Weizsacker-Williams (WW) distribution [11, 44-47]. 
This approach has been widely used with much success in the 
HERA physics program [11]. 

In the context of our NLO calculation for iN hX it is 
therefore interesting to investigate whether also in this case 
the contributions by almost real photons dominate and the 
NLO corrections may be well approximated by a Weizsacker- 
Williams type distribution. Since it is much easier to com¬ 
pute the latter contribution than the full NLO correction, this 
would mean that one could also obtain approximate NLO re¬ 
sults for the transversely polarized cross section within the 
twist-3 framework by simply considering real photons. Given 
the complexity of a full NLO calculation for the twist-3 case, 
this would be a tremendous advantage. We note that the con¬ 
tributions to the spin-dependent cross sections for iN —> jet Y 
for real photons were discussed in [6], including the twist-3 


contributions for the single-transverse spin case. Actual LO 
calculations for the twist-2 longitudinal spin-dependent cross 
section were presented in Ref. [48] for quasi-real photons. We 
will closely examine the contributions by quasi-real photons 
also in our paper. Their relevance will of course also depend 
on the lepton species that is used, because the lepton mass 
leads to a lower limit on the virtuality of the photon. 

Our paper is structured as follows. In Sec. II we present our 
NLO calculations for the partonic cross sections for iN'^ —> 
hX and iN —> jetX. We also discuss in some detail the 
Weizsacker-Williams contribution and how the calculation 
can be done keeping a finite lepton mass. Section III presents 
numerical predictions for the NLO cross section to be ex¬ 
pected at various fixed-target experiments and at a future Elec¬ 
tron Ion Collider (EIC). Einally, we summarize our results in 
Sec. IV. 


II. NLO CALCULATION 
A. General framework 

In this section we present our derivation of the analytical 
NLO results for the processes iN —» hX and iN —» jet A. 
The transverse momentum of the produced hadron or jet sets 
a hard scale, so that perturbative methods may be used for 
treating the cross sections. We first consider i{l) -t- N{P) 
h(Ph)+X, where we have introduced our notation for the four- 
momenta. It is useful to introduce the Mandelstam variables 
asS - (P + l)^, T - (P-Ph)^ and U - {l-Ph)^. Eurthermore, 
we label the energy of the detected hadron as £/, and its three- 
momentum by Ph ■ 

In collinear leading-twist perturbative QCD the hadronic 
cross section is approximated by convolutions of hard partonic 
scattering cross sections and parton distribution/fragmentation 
functions. The momenta of the incoming parton, and of 
the fragmenting parton, fP, which appear in the calculation of 
the partonic cross sections, are approximated as ^ xP^ and 

- P^ jz, respectively. It is then convenient to work with the 
partonic Mandelstam variables 

s-{k + l)^-xS, t-{k-p)^--T, M = (Z - p)^ = — . (1) 

z z 

The general form of the factorized cross section for the inclu¬ 
sive hadron production process then is 



X D'’^^(z,p) d-‘^^(s,t,u,p) , (2) 

where f‘^^(x, p) is the parton distribution function (PDF) for 
the incoming parton i in the nucleon N and D''^^\z,p) the 
corresponding fragmentation function (FF) for parton / frag¬ 
menting into hadron h, both evaluated at a factorization scale 
p. We choose the factorization scales to be the same for the 
initial and the final state, and also equal to the renormaliza¬ 
tion scale. In Eq. (2), is the partonic cross section for 
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FIG. 1. LO diagram for lepton-quark scattering. 


the lepton-parton scattering process, { + i f + x, with x 
an unobserved partonic final state. The sum in Eq. (2) runs 
over the different species of partons, quarks, gluons and an¬ 
tiquarks. We note that the expression in Eq. (2) holds up to 
corrections that are suppressed by inverse powers of the pro¬ 
duced hadron’s transverse momentum Pi,j_. 

The partonic cross sections in Eq. (2) can be calcu¬ 
lated in QCD perturbation theory. One may write their expan¬ 
sion in the strong coupling as 


cr 


= cr, 


Hf 

LO 


+ — cr, 


Hf 

NLO 


H- 0{a]). 


(3) 


At lowest order (LO) only the tree-level process {q q( 
shown in Eig. 1 contributes. The calculation of its cross sec¬ 
tion is straightforward. One finds 


cr, 


LO 


- 2a e 




6{s + t + u) , 


(4) 


where Ogm is the fine structure constant and is the quark’s 
fractional charge. 

At NLO, both virtual (Eig. 2) and real-emission 

diagrams (Eigs. 3a-3c) contribute. We will address these in 
turn in the following subsections. One can see from Eigs. 3b 
and 3c that beyond LO there are also new contributions where 
a gluon fragments or where an initial gluon enters the hard 
scattering process. 

As is well known, all types of NLO contributions develop 
singularities at intermediate stages of the calculations, which 
we make manifest by using dimensional regularization with 
D = 4 - 2e space-time dimensions. The subsequent treatment 
of the singularities is standard in pQCD calculations. The only 
non-standard feature arises for the incoming lepton. If we as¬ 
sume for the moment that we have an incoming quark instead 
of a lepton in the diagrams in Eigs. 3a and an exchanged gluon 
instead of a photon, then the diagram would make an NLO 



FIG. 2. Virtual diagrams at NLO. Self energy diagrams (right and 
middle graph) contribute in Feynman gauge. 


contribution to, say, pp —> hX. Being treated as massless, the 
initial quark would produce a singularity when it radiates the 
gluon collinearly. As is well understood, this singularity may 
be absorbed (“factorized”) into the proton’s quark PDF, ex¬ 
actly in the same way as for the incoming quark at the bottom 
of the diagram. In case of an incoming lepton, on the other 
hand, the lepton’s mass ensures that no collinear singularity 
arises when the lepton radiates a collinear photon that subse¬ 
quently participates in the hard scattering. In fact, keeping 
the lepton mass m(, the cross section will develop a logarith¬ 
mic term of the form aem log(A/m^), where A represents a 
hard scale of the problem, and in the limit m( —> 0 this log¬ 
arithm precisely produces the required collinear singularity. 
In principle we should therefore perform the NLO calculation 
keeping the lepton mass finite. This is technically very cum¬ 
bersome, and in fact not needed. We can adopt two different, 
and equivalent, approaches instead: In the first approach we 
neglect the lepton’s mass and regularize the ensuing collinear 
pole in dimensional regularization. The pole is then subtracted 
(for example, in the MS scheme) and absorbed into a “parton” 
distribution function for photons in a lepton. This distribution 
may be evaluated perturbatively in first-order QED, giving rise 
essentially to the well-known “Weizsacker-Williams” distri¬ 
bution. This approach may in principle be extended to higher 
order in QED. In the second approach, we calculate the cross 
section for a massive lepton, keeping however only the leading 
terms in nif which are of the form aem[iog(A/nif) + constant]. 
This is justified by the fact that all terms beyond this approx¬ 
imation are suppressed as powers of m( over the hard scale 
and hence numerically tiny. We note that although the loga¬ 
rithm can become large (as me is small compared to typical 
QCD hard scales), the smallness of aem will usually make the 
term a log(A/OTf) small enough to be regarded as a perturba¬ 
tive correction. We will present our main calculation for the 
case of massless leptons and comment on the use of a finite 
lepton mass in the calculation later. 

It is convenient to rewrite the x- and z-integrals in Eq. (2) 
in terms of new variables v - 1 + tjs and w - -u/{s + f). 
Using (1), we have 

l-vU -T 

X = -^ = n -^ 

vw 1 (1 - v)o 


and Eq. (2) becomes 




r'+f 

dv r' dw 


^ d^Ph 

\S^ 


v(l - v) 

V T 



X H‘Hv,w) 


W, jl) , 

(6) 

where we have defined 




H'f{v,w) = 

X z? 

1-y U -T 

(7) 


vw 


For ease of notation, we have kept the symbol also for 
the cross section when expressed in terms of the new vari¬ 
ables. We note that the invariant mass of the unobserved re¬ 
coiling partonic final state is given hy s + t + u - sv(l - w). 
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(a) 


(b) 


(c) 


FIG. 3. NLO real-emission diagrams. There are three partonic channels at NLO: (a) ^ (b) g’ —> g, (c) g —> g. 


The function -H f H- m) oc 5( 1 - w) in the LO cross section (4) 
expresses the fact that at LO the recoil consists of a single 
parton. 


in (10) requires special care in this limit since the denominator 
would lead to a non-integrable behavior for e = 0. We deal 
with this limit by means of the expansion 


B. Virtual contributions at NLO 


At the NLO level, the virtual contributions shown in Fig. 2 
contribute through their interference with the Born diagram. 
The virtual contributions thus have Born kinematics and are 
proportional to 5(1 - w). Since we are only interested in QCD 
virtual corrections, only the quark line is affected, and we may 
adopt the result directly from the corresponding calculation in 
Ref [49] for the basic photon-quark scattering diagrams in 
DIS. This gives 

„^ Cfa,ip) r(i -g)^r(i -Hs) 

O'NLO.vir - 27T r(l - 2s) 


X 



where 




_ 

LO,f 


= 2a; 


2 ^ 
e — 

em £7 
^ SV 


/ 1 H- 
1(1 - V )2 



5(1 - w) . 


( 8 ) 

(9) 


(1 - 


1 

2s 


5(1 -w) + 


1 

(1 - w)+ 


- 2s 


ln(l -w)\ 
l-w 


+ 0{S^) , 


( 11 ) 


where the plus distribution is defined in the usual way by 

f dw f(w) lg(w)]_^ ^ f dw [f(w) - f(l)] g(w) . (12) 

Jo Jo 

This expansion makes the singularities in 1 js explicit. When 
combined with the pole terms in the term oc 5(1 - w) 

in (11) leads to a double pole term that cancels against the 
double pole in the virtual correction in Eq. (8). This well- 
known behavior reflects the cancelation of infrared singulari¬ 
ties in partonic observables. The channels q ^ g and g ^ q 
in Figs. 3b and 3c are infrared finite at NLO. 


is the Bom cross section computed in 4 - 2e dimensions. Fur¬ 
thermore, Cf - (NI — \)l2Nc, with Nc the number of colors. 


C. Real-emission corrections at NLO 


D. Collinear subtraction for parton distribution functions and 
fragmentation functions 


The real diagrams have 2^3 topology. To obtain the 
desired contribution to an inclusive-parton cross section we 
need to integrate over the phase space of the lepton and the 
“unobserved” parton in the final state. This can be done in 
4 - 2e dimensions using the standard techniques available in 
the literature [50-52]. 

After phase space integration, the result for the real- 
emission contribution for the q ^ q channel takes the form 


O-Nrareal //,£) + 


W,IJ.,s) 

(1 _ w)l+2fi 


( 10 ) 


where both functions and d"^^^ carry a 1 /e-pole, but are 
well-behaved in the limit w ^ 1. Obviously, the second term 


After the cancelation of infrared singularities between real 
and virtual contribution, the partonic cross sections still ex¬ 
hibit single poles that reflect collinear singularities arising 
when an “observed” parton (either the incoming one, or the 
one that fragments) becomes collinear with the unobserved 
parton. The factorization theorem states that these poles may 
be absorbed into the parton distribution functions or into the 
fragmentation functions. This procedure may be formulated 
in terms of renormalized parton densities and fragmentation 
functions (see, e.g.. Ref. [53]). In fact, naive definitions of 
“bare” parton densities and fragmentation functions contain 
ultraviolet singularities that can be dealt with as well by using 
dimensional regularization. At NLO, the corresponding ultra¬ 
violet 1/e-poles that appear can be removed in the MS scheme 
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by introducing “renormalized” functions in the form 


/ble = /ref ix,n)+ ^[Pqq® ffjn) 


S e 
In e 


[Pq,®f!l^)(x,ii) + 0{a^;), (13) 


^ {Pqq ® £>ren) (z,^) 


27r s 


S £ 


27!' S 

where we have the usual splitting functions 


{Pgq®D>ii^)(z,^i) + 0{al), (14) 


Pqqiy) = Cf 


T ^^ + ^rJd-y) 

(i-y)+ 2 


Pqsiy) = + ’ 

l+(l-y)2 


Pgq(y) — Cf 


y 


(15) 

(16) 
(17) 


(with r« = 1 /2), and where the ”®”-symbol indicates the con¬ 
volution 


=r 


dy 


{P®f)ix)= -P{y)f\- 


y 


(18) 


The constant = (47r)®/r(l -s) in (13) and (14) corresponds 
to the usual MS scheme. Inserting the bare distributions into 
the LO expression for the hadronic cross section, we obtain 
additional 0{a,,al^) contributions. These precisely cancel the 
collinear poles associated with the observed partons in the 
NLO partonic cross sections, for all three channels. 

Even after this procedure, one type of collinear singularity 
remains. It is generated by a momentum configuration where 
the exchanged photon is collinear to the incoming lepton. As 
discussed at the beginning of this section, the presence of this 
singularity is an artifact of neglecting the lepton’s mass. In 
the following two subsections we discuss our treatment of this 
issue. 


E. Weizsacker-Williams contribution 

One approach for dealing with the collinear lepton singu¬ 
larity is to introduce bare and renormalized QED parton dis¬ 
tributions for the lepton, much in analogy with the procedure 
that we discussed in the previous section for the nucleon’s par- 
ton distributions. The only differences are that for leptons the 
partons are the lepton itself and the photon, and that we can 
safely compute their distributions in QED perturbation theory. 
To lowest order in QED, we have just f^^^(y) - b(l -y), corre¬ 
sponding to the Born contribution in Eig. 1. The hard process 
involving an incoming lepton will always require two electro¬ 
magnetic interactions and hence be of order as seen ex¬ 
plicitly in Eq. (4). This is different for a hard process with an 
incoming photon such as yq qg, which is of order aemO's- 
This implies that at NLO in QCD (at order there will 



FIG. 4. General Weizsacker-Williams contribution at NLO. The 
quasi-real photon entering the hard scattering part is treated as a par- 
ton in the lepton. 


be contributions generated by the photon acting as a parton 
of the lepton and participating in the hard process. A generic 
picture for such types of contributions, known as Weizsacker- 
Williams contributions, is shown in Eig. 4. In essence, the 
lepton merely serves as a source of real photons for the con¬ 
tributions shown in the figure. Like its nucleon counterpart, 
the corresponding photon-in-lepton distribution P^^iy) will 
require renormalization. Eollowing (13) we may write 

/toeCy- p) ^ flJn(y,P)+‘^^ [Py( ® /re(f) O'- • (19) 

where Pyt - PgqlCf and the ellipses denote a term involving 
a photon-to-photon splitting that makes contributions beyond 
the order in aem we consider here. Within the same reasoning, 
we can set ////(y) = <5(1 - y) in (19). ^ 

The bare photon-in-lepton distribution in Eq. (19) can 
be defined analogously to the gluon distribution in a nucleon 
in terms of the matrix element (see also [6]) 

Q''^(y) = npn^ f ^ (eiFZ(0) U[0; An] FZ(MP, 

J-ac 27ry 

= (20) 

In this definition « is a light-cone vector conjugate to the lep¬ 
ton momentum /, with n^ - 0 and I ■ n - 1. Eurthermore, 
^em = - d’'A^‘ is the electromagnetic held-strength ten¬ 

sor, and we have inserted a (straight) Wilson line U[Q',An] 
that ensures the electromagnetic gauge invariance of the ma¬ 
trix element. The transverse projector in (20) is given as 
gZ = g^‘'' - Pn'' - 

Since the matrix element in (20) contains electromagnetic 
fields and elementary leptons in the in- and out-states we can 
compute it to LO in QED. In this calculation we keep a non¬ 
vanishing lepton mass m{ in order to obtain an infrared-finite 
result. To order Oia^ra) we hnd. 




[iH-ln 

f 1 

- 1 

e 

[y^mj) 



■ </(«em) > 

( 21 ) 


where, as before, S e = (47r)®/r(l - e). In close analogy to 
parton distributions of the nucleon we can perform an MS- 
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renormalization of the distribution and obtain. 


' bare ' 


Q'em n / N 

17 


^em 

In 

Pyc(y) ^ 
e 


( 2 \ ^ 

In 

, -1 


[y^mj) 


+ 0(0. (22) 


This renormalized distribution is closely related to the ‘clas¬ 
sic’ Weizsacker-Williams distribution [44^7]. The logarithm 
in (22) may be derived from an integration over the pho¬ 
ton’s virtuality -q^ (where q is the photon momentum). For 
the standard Weizsacker-Williams distribution one performs 
this integration from the lower kinematic limit - y) 

to an upper limit fixed by the experimental condi¬ 

tion imposed on the scattered lepton. This gives rise to a 
term ^Pjciy) ln(2^^(l -y)l(y^m])) in the photon spectrum, 
which can be recovered by an appropriate choice of the scale 
p in (22). 

For the contribution related to the photon virtuality is 
then neglected everywhere else in the hard scattering. One 
thus considers scattering diagrams with a real incoming pho¬ 
ton. We thus write the generic factorized cross section for the 
contribution as 



X f^Jn\y,n) (23) 

with the cross sections describing the scattering yi 

fx of the photon off parton i in the nucleon (to be given be¬ 
low). At 0(as) we encounter three channels with an incoming 
photon: yq q{g), yq —> g{q), and yg q{q) (the partons in 

parentheses are not observed). The relevant diagrams are as 
those shown in Figs. 3a - 3c, but with the lepton lines removed 
and the virtual photon replaced by a real photon. Being 2 —» 2 
diagrams, their calculation is straightforward. Inserting now 
the bare WW distribution we generate precisely the pole terms 
required to cancel the lepton collinear divergences discussed 
at the end of Sec. IID. This happens in the same way for all 
partonic channels. We note that the dependence on the scale fi 
associated with the lepton also disappears. This has to be the 
case, since for a finite lepton mass there would never be any 
lepton collinear divergences in the first place. 


F. Calculation with nif t 0 

As we noted earlier, the presence of collinear singularities 
associated with lepton-photon splitting is really an artifact of 
neglecting the lepton’s mass. In principle we should there¬ 
fore perform a full calculation in which the lepton’s mass is 
kept finite. This is trivial for the virtual diagrams, since the 
QCD corrections do not affect the lepton line. However, in¬ 
clusion of a lepton mass considerably complicates the phase 
space integrations for the real diagram. Nevertheless, it is pos¬ 
sible to compute the relevant integrals using the results given 
in Ref. [51]. One may then expand the result in powers of the 


lepton mass and neglect terms suppressed by powers of Oimi). 
In this way, the “would-be” collinear singularity is regularized 
by the lepton mass and shows up as a term ~ In(m^). Terms 
independent of nii are also kept. All other parts of the calcula¬ 
tion proceed as before, and the partonic cross section thus has 
the structure 

^ ln(m^/s)-H 

w,ii) + 0(m^ ln(m^)). (24) 

for each channel. 

We have checked explicitly for all three channels that our 
two approaches for treating the initial lepton are equivalent: 
The full result obtained using the WW contribution in the 
previous subsection agrees with that for m{ + 0, as long 
as we only keep the leading terms as discussed in Eq. (24). 
The equivalence of the two approaches serves as an important 
check of our calculation and also explicitly demonstrates the 
universality of the WW-distribution. 


G. Final results for single-inclusive hadron production 


We now present our final results for the full partonic cross 
sections in analytic form. Combining the cross section (6) 
for massless leptons with the Weizsacker-Williams contribu¬ 
tion (23), we may write the full NLO cross section as 


Eh- 


cPcr^^ 


■hx 


d^Ph 



dv 1 

r' dw 

'52 jp J 

^ v(l-v)J 

'l-vU ^2 

V T 


X H‘d(v,w) 


'^LO iv) + 


fy/f ( l-v A N 

Jren \ ^ ^TlO '^1 


(25) 


where Hd(v, w) has been defined in Eq. (7). The LO contri¬ 
bution, present only for the channel q ^ q with an incoming 
quark that also fragments, was already given in (4). Eor the 
NLO term in this channel we find 




q.2 


1 


(1 - w)+ 


Apd(l-w)+Ap 


ln(l - w) 
1 - w 


v(l - v(l - w)) 


sv 


(1 - v)w 


+ ‘'ln(y(l - w))-hC: 

1 - V 


pin 


1 - vw 


H- Cpin 


(1 - vw)(l - v(l - tv)) 


+ cr^n ^ +c: 


(26) 
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where the coefficients are functions of v and 

w and may be found in the Appendix. The channels q ^ g 
and g q have simpler expressions: 




allele F 


+ C^'^ln 


+ Cf>^ln 


Cp^ln(l - v(l - w)) 
1 -V 

(1 - vw)(l - v(l - w)) 
v(l - w)s 




^ *-4 


(27) 


^- 4 , , ^LelTR 


svw 


+ Cf^^ln 


C^^ln 


v(l - w)i 


(1 - v)w 
1 - vw 




(28) 


The coefficients and are again given in the Ap¬ 

pendix. 

We finally list the partonic cross sections for the 
Weizsacker-Williams contributions: 


LO 


LO 


cr\ 


LO 


(v, w) 

(v,w) 

(v, tv) 


2,71 C 1 -t V^tV^ 
i( 1 - V) VW 

InCpa^mel 1 h- (1 -vwf 
i(l - v) 1 - vw 

2n TRUf^rne^ v^w^ + (1 ~ vw)^ 
i(l - v) vw(l - vw) 


(29) 


H. Single-inclusive jet production 


Having computed the inclusive hadron production cross 
section at NLO the extension to single inclusive jet produc¬ 
tion is straightforward. The cross section for {N —»jetX may 
be written as 




d^Pj 


' ^ S+T 


X (v = 1 -I- w,;u;, (30) 

where Ej and Pj are the energy and three-momentum of the jet 
and the hadronic Mandelstam variables are defined as before, 
now in terms of the jet momentum. The form of this expres¬ 
sion follows from (6) by setting the fragmentation functions 
to (5(1 - z). Of course, beyond LO, the partonic cross sec¬ 
tions for jet production differ from the ones for single- 

inclusive hadron production. This is evident from the fact that 
the latter are computed as “inclusive-parton” cross sections 
which, as we saw in subsection IID, require collinear 
subtraction. This is in contrast to a jet cross section which is 
by itself infrared-safe, as far as the final state is concerned. 
Instead, it depends on the algorithm adopted to define the jet. 


as we have indicated by the dependence on a generic jet (size) 
parameter R in (30). 

As was discussed in Refs. [54-56], even at NLO one may 
still go rather straightforwardly from the single-inclusive par- 
ton cross sections to the for any infrared-safe jet 

algorithm. The key is to properly account for the fact that 
at NLO two partons can fall into the same jet, so that the jet 
needs to be constructed from both. In fact, assuming the jet to 
be relatively narrow, one can determine the relation between 
and analytically [54]. This “Narrow Jet Approx¬ 
imation (NJA)” formally corresponds to the limit B —> 0 but 
turns out to be accurate even at values B ~ 0.4 - 0.7 rele¬ 
vant for experiment. We follow this approach in this work. 
In the NJA, the structure of the NLO jet cross section is of 
the form ,^Iog(B) + S; corrections to this are of 0{R^) and 
are neglected. We note that to the order we consider 

in this paper, the Weizsacker-Williams terms only contribute 
to the B-independent piece S. This is because for almost real 
exchanged photons it is at this order not possible to have two 
coalescing partons in the final state. 


III. NUMERICAL RESULTS 


We now present phenomenological results for the NLO 
single-inclusive pion production cross section in lepton- 
proton scattering. As mentioned before, data on the transverse 
single-spin asymmetry for this reaction have been released 
by HERMES [2] and the Jefferson Lab Hall A Collabora¬ 
tion [4]. Unfortunately, corresponding cross sections were not 
presented, and we will therefore provide predictions for these. 
Eurthermore, we will also present predictions for COMPASS 
at CERN, for a future Electron-Ion-Collider (EIC), and for 
experiments at Jefferson Lab after the CEBAE upgrade to 12 
GeV beam energy. Einally, at the end of this section we show 
some phenomenological results for the inclusive production 
of jets at the EIC. 

As we saw in the previous subsections (see Eq. (25)), our 
NLO result can be formulated in such a way that it contains 
contributions involving the photon-in-lepton distribution 
and LO photon-parton cross sections. These represent the 
contributions by quasi-real photons to the cross section. An 
interesting question is whether this part of the cross section 
dominates the NLO corrections, at least for a suitable choice 
of the scale in (22). We recall that the logarithm in (22) 
may be obtained by an integration over the photon’s virtuality 
where only the Ijq^ propagator is kept for the photon, while 
q^ is neglected everywhere else in the hard scattering. We now 
consider the cross section 


Eh 


d^Ph 


l~U\'SA dv r' dw 

v(l - v) Jtzz E w^ 

ij T+u ' V T 


which essentially corresponds to the full NLO one in (25), 
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FIG. 5. Cross section for {p —> n^X at HERMES, (a) as function of Xf for 1 GeV < Pi,i_ < 2.2 GeV, and (b) as function of for 
0.3 < Xf < 0.55. The dashes line gives the LO prediction and the solid line the NLO one. The dotted and dot-dashed lines show the 
approximation (31) of the NLO cross section, using pq = Pi,f and pq = VS/2, respectively. 


but with the terms dropped. In other words, we use 

the LO term and add the Weizsacker-Willia ms co ntribution. 
For the latter, we choose the upper limit on yj-q^ in the pho¬ 
ton spectrum as a large scale in the problem, pQ ~ or 
even /tq ~ VS12. This constitutes an attempt to obtain an 
approximation to the full NLO correction by assuming that 
the 1 -behavior of the hard cross sections is valid over most 
of the kinematical regime. In our studies we examine in this 
way the importance of the Weizsacker-Williams contribution. 
As discussed in the Introduction, if the contribution plays a 
dominant role for the NLO corrections, this opens the door 
to approximate NLO calculations also for the spin-dependent 
case. 

For all our calculations we use the CTEQ6.6M [57] set of 
parton distribution functions and the fragmentation functions 
of [58]. 


1. HERMES 


Figures 5 a and 5b present our results for -production at 
HERMES at yfS — 7.25 GeV. We fix the renormalization and 
factorization scales at ju = Phi_. In order to match the conven¬ 
tions used mpp^ —> hX, HERMES presents the spin asymme¬ 
try results in terms of the hadron’s transverse momentum 
and Eeynman’s xp - 2F~/ xfS, where is the z-component 
of the hadron momentum in the center-of-mass frame of the 
collision, and where positive is counted in the direction of 
the lepton beam. We have 


dxf dPhi. 


2nPhi. 

+ 4 


Eh 


fig-ep-^T^X 

d^Ph 


(32) 


where xp - 2Ph±l sfS. The hadronic Mandelstam variables 
read 

T = -|(74 + 4+'^f) ’ 

U = +xl- xf^ . (33) 

Figure 5a shows the cross section as a function of xp, inte¬ 
grated over 1 GeV < F/,x < 2.2 GeV. This is the only Phi. bin 
used in Ref. [2] with Phi > 1 GeV. In Fig. 5b we examine 
the Phi dependence of the cross section for 0.3 < xp < 0.55. 
In both cases we find large NLO corrections; the NLO cross 
section is almost twice as large as the LO one. As discussed 
above, we also examine in how far the Weizsacker-Williams 
contribution drives the NLO corrections, using Eq. (31) with 
Fo = Phi (dotted) and j2 (dot-dashed). As one can 

see from the figures, the Weizsacker-Williams contribution 
does lead to an increase over LO, but provides only about 50% 
to 70% of the NLO correction. This is likely to be attributed to 
the fact that the overall c.m.s. energy is rather low. The result 
with pq = sfS /2 provides a slightly better description of the 
full NLO, although the differences are minor. We note that the 
WW approximation appears to work better for smaller trans¬ 
verse hadron momenta Phi and for larger xp. The latter fea¬ 
ture perhaps is at first sight surprising since positive xp of the 
hadron imply on average backward scattering of the lepton, 
whereas the WW approximation should work better if the lep¬ 
ton is scattered in the forward direction. One can roughly un¬ 
derstand this “shift” of the WW approximation towards posi¬ 
tive Xp from the fact that |r| » |t/| for .x/r ^ 1 in Eq. (33). 
Since the dominant real-photon process yq q(g) in (29) has 
a 1 / iM-behavior in contrast to the 1 /f^-behaviour of the LO 
process, the WW aproximation favors the region xp > 0. The 
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(a) (b) 

FIG. 6. Same as Figs. 5a, 5b, but for ^ ^He scattering at beam energy 12 GeV after the CEBAF upgrade at Jefferson Lab. For the cross section 
as a function of Xp we have used a fixed = 1.5 GeV, while for the Pjx dependence we have integrated over -0.4 < Xp < 0.4. 


full NLO partonic cross section inherits the 1/f^-behaviour of 
the LO one, so that the Weizsacker-Williams contribution can 
approximate it well only for xp > 0. 


2. Scattering with the 12 GeV beam at the Jefferson Lab 


Our NLO predictions for the cross section for ^ ^He ^ 
in 12 GeV scattering at the Jefferson Lab are shown in Figs. 6a 
and 6b. For the xp distribution on the left we have assumed a 
fixed transverse momentum Ph± = 1.5 GeV. On the right we 
show the Pi,j_ dependence of the cross section in the region 
-0.4 < xp < 0.4. Again, the renormalization scale is fixed 
to the transverse hadron momentum, fi - Phj_. Note that the 
rather modest c.m.s. energy available limits the possible size 
of P/,j_ severely. For collisions using the present 6 GeV beam 
only transverse momenta outside the hard-scattering regime 
are possible, which is the reason why we cannot present any 
results for this case. 

We again observe in Figs. 6a, 6b that the NLO correc¬ 
tions are very large. The Weizsacker-Williams contribution 
is clearly insufficient to match the NLO result here. 


3. COMPASS 


forward direction of the incident muon. We have 

- = 2nPf,j_Ef, - 5 - , 

dtjdPh^ d^Phi. 

where the hadronic Mandelstam variables read 

T = -P/,^V5e+^ 

U = -Ph^sfSs-’’. 


(34) 


(35) 


The COMPASS spectrometer roughly covers the region 
-0.1 < T] < 2.38. From the t] dependence shown in Fig. 7a 
for a fixed transverse momentum Ph± - 2 GeV we observe 
that the NLO corrections are significant but not as large as for 
HERMES and JLab. They amount to an increase over LO of 
roughly 30-40%. Strikingly, the Weizsacker-Williams contri¬ 
bution is very small here, even for the choice po = sfS /2. This 
may be understood from the fact that the muon mass is about 
200 times larger than the electron mass, resulting in a much 
smaller logarithm in the expression (22) for the photon spec¬ 
trum, which then is largely cancelled by the non-logarithmic 
term. 

Eor the Ph± spectrum shown in Fig. 7b we also show the 
results for a different choice of the factorization and renor¬ 
malization scales, p = 2Phj_. As one can see, the scale depen¬ 
dence decreases somewhat when going from LO to NLO but 
remains fairly sizable. 


The results of our NLO analysis for COMPASS kinemat¬ 
ics are shown in Figs. 7a and 7b. COMPASS uses a muon 
beam with energy 160 GeV, resulting in sfS - 17.4 GeV. 
Following the choice made by COMPASS, we use here the 
c.m.s. pseudorapidity t] of the produced hadron rather than 
its Feynman-x/r. Pseudorapidity is counted as positive in the 


4. Electron-Ion Collider 

We finally also discuss the cross section for single-inclusive 
pion production in electron-proton collisions at a proposed fu¬ 
ture EIC [59] with yfS - 100 GeV. Thanks to the higher en¬ 
ergy of an EIC it will become possible to probe much larger 
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(a) (b) 

FIG. 7. Cross section for /rp —» tt^X at COMPASS, (a) as function of hadron pseudorapidity for fixed P;,j^ = 2 GeV, and (b) as function of P/,j^ 
for -0.1 < Tj < 2.38. As before, the solid lines give the NLO results and the dashed lines the LO ones. The dotted and dot-dashed lines show 
the approximation (31) of the NLO cross section, using fio = P|,^^ and fiQ = yfS 12, respectively. In (b) we also present the LO and NLO results 
for the scale = 2Pi,±. 


transverse hadron momenta, where pQCD is expected to work 
better. Fig. 8 a shows the rj dependence of the cross section for 
a fixed transverse momentum P/,^ = 10 GeV. Again we count 
positive 77 in the forward direction of the incoming lepton. The 
P/,x dependence of the cross section is shown in Fig. 8 b, in¬ 
tegrated over It/I < 2. The renormalization scale has again 
been fixed to the transverse hadron momentum, /j. - Pf,±. As 
for COMPASS we found that the scale dependence slighty de¬ 
creases for EIC kinematics when going from LO to NLO but 
remains relatively large. 

We again find sizable NLO corrections. Overall, the 
Weizsacker-Williams approximation works much better here 
than in the fixed-target regime. It describes the NLO cross sec¬ 
tion especially well when the hadron is produced in the elec¬ 
tron forward direction. At mid-rapidity and negative rapidity 
the approximation tends to fall short of the full NLO result. 
Lrom Lig. 8 b we observe that that the Weizsacker-Williams 
also works better for smaller 


5. Jet production at an EIC 

Given the high energy of an LIC, also jet observables will 
be of much interest there [ 6 ]. Lor example, combined analy¬ 
sis of data for the transverse-spin asymmetries for ep'' —» hX 
and ep^ jetV from a future LIC should allow for a clean 
separation of twist-3 parton correlations in the nucleon and in 
fragmentation. We therefore close this section by presenting 
predictions for the cross section for single-inclusive jet pro¬ 
duction, ep —> jet V. Here we use the NJA formalism outlined 
in Sec. IIH to convert the single-hadron cross section into a jet 
one. We adopt the anti-k, jet algorithm of [60]. In Lig. 9a we 


present the dependence of the cross section on the jet pseudo¬ 
rapidity 7 /j for a fixed transverse jet momentum of Pj^ = 10 
GeV. We find once again that NLO contributions are large. We 
also observe that, compared to the case of hadron production 
considered in Lig. 8 a, the NLO cross section is much more 
peaked in the forward electron region. The reason is that at 
large positive pseudo-rapidity |r| » |t/| in Lq. (35). Since 
the minimal value for the incoming parton’s momentum frac¬ 
tion is Xjnin = -UI{S + T) in (30) rather small values of x are 
probed at large pseudo-rapidity where in turn the nucleon’s 
parton distributions are large. On the other hand the fragmen¬ 
tation process suppresses the forward and backward regions in 
hadron production due to the large z-values probed, whereas 
in jet production the forward electron region is enhanced due 
to the absence of fragmentation. 

In the figures 9a and 9b, we show results for two different 
jet size parameters, R = 0.7 and R - 0.2. Dependence on R 
first occurs at NLO. As discussed at the end of Sec. IIH, the 
first-order Weizsacker-Williams contribution does not depend 
on R. It hence cannot give an accurate approximation of NLO 
in general. As the figure shows, the WW result happens to be 
rather close to the result for R - 0.7; this agreement, however, 
is essentially fortuitous. 


IV. CONCLUSIONS AND OUTLOOK 

We have performed next-to-leading order calculations of 
the partonic cross sections for the processes {N hX and 
JN jet V, for which the scattered lepton in the final state is 
not detected. We have derived our results for a finite lepton 
mass, neglecting terms that are suppressed as powers of the 
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(a) (b) 

FIG. 8 . Cross section for ep —> n^X at an EIC with yfS = 100 GeV, (a) as function of p at fixed = 10 GeV, (b) as function of 
integrated over I 77 I < 2. The lines are as in the previous figures. 


mass over a hard scale. The results have been obtained in two 
ways. We have first set the mass to zero. We have regularized 
the ensuing collinear singularity in dimensional regularization 
and then subtracted it by introducing a Weizsacker-Williams 
type photon distribution in the lepton. The latter can be com¬ 
puted in QED perturbation theory and effectively reinstates 
the leading lepton mass dependence, which is logarithmic plus 
constant. In the second approach, we have kept the lepton 
mass in the calculation directly, expanding all phase space 
integrals in such a way that the leading mass dependence is 
obtained. Both approaches give the same result. 

We have presented phenomenological NLO predictions for 
various experimental setups, from fixed-target experiments 
(HERMES, JLab, COMPASS) to collider experiments at an 
EIC. We have found that the NLO corrections are large. We 
note that in the fixed target regime the bulk of the corrections 
comes from the plus distribution terms in Eq. (26), especially 
at negative ;ic/r or rapidity. As is well known, the distributions 
are associated with the emission of soft gluons. Since they 
recur with increasing power at every higher order of perturba¬ 
tion theory, it may be worthwhile for future work to address 
their resummation to all orders, similarly to what was done for 
the photoproduction case £N —> fhX in [61]. 

The rather large size of the corrections that we find sug¬ 
gests that also the cross section with transverse polarization 
of the initial nucleon may be subject to large NLO correc¬ 
tions. This would likely have ramifications for analyses of 
spin asymmetry data for hX in terms of twist-3 parton 

correlation functions. As full NLO calculations for transverse 
single-spin observables are difficult, we have also investigated 
in how far it is possible to match our full NLO result for the 
spin-averaged cross section by adding just the Weizsacker- 
Williams contribution to the LO one. We have found that this 
simplified approach does not appear to work well quantita¬ 


tively. In other words, the NLO corrections do not appear 
to be dominated by quasi-real photons. Nonetheless, in or¬ 
der to obtain a first estimate of higher-order effects for the 
transverse-spin asymmetry, it may be worthwhile to use the 
Weizsacker-Williams contribution for the case of transversely 
polarized nucleons, which is much simpler to do than the full 
NLO calculation and was already discussed in Ref. [6]. 

We again emphasize that our results suggest that contribu¬ 
tions by quasi-real photons to the cross sections for the single- 
inclusive processes {N —> hX and (N —> jetX are not the 
dominant contributions, at least for large transverse hadron 
momenta > 1 GeV. In other words, an experimental 
setup where the final state lepton is not observed in lepton- 
nucleon collisions does not automatically imply that one mea¬ 
sures an (approximated) quasi-real photoproduction process. 
However, although quasi-real photons do not dominate, they 
typically do play a non-negligible role for the NLO correc¬ 
tions. As is well known, high-energy real photons may also 
exhibit their own partonic structure, in which case they are 
referred to as “resolved” photons (see Ref. [11]). The cor¬ 
responding resolved-photon contributions are formally of the 
same order as the Weizsacker-Williams contribution we have 
considered here. They are typically suppressed in the fixed- 
target regime. It may be interesting to address this contribu¬ 
tion in future work, also in order to study its impact on the 
transverse single-spin asymmetries. The concept of “virtual 
photon structure” may also prove useful in this context (see, 
for example. Ref. [62]). 
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(a) (b) 

FIG. 9. Cross section for single-inclusive jet production at the EIC, (a) as a function of pseudorapidity r/j at a fixed transverse jet momentum 
Pjx = 10 GeV, and (b) as function of integrated over |77 j| < 2. We have used the NJA [54, 55] and the anti-fe, jet algorithm [60]. The solid 
and dotted lines show NLO prediction for two different values of the jet size parameter, R = 0.7 and R = 0.2, respectively. The dashed lines 
present the LO results, and the dotted ones the result for the approximation (31) of the NLO cross section, using yUo = P!± and (on the left) also 
Ro= yfS/l. 
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Appendix A: NLO coefficients 

Here we present the NLO coefficients in Eqs. (26), (27), 
(28) for the different channefs; 
a. q ^ q channel: 


A‘1-^9 

'^0 








^(3-b 2fn(v))fn( '^^^ -H fn^(v) - 8 I, 


= 8 w 


= 4w 


(l-v)2 

1 -H v 2 




(1-V)2’ 

1 - v(l - w) H- v2(l - w(l - w)) 


(1 - V )2 


2 w 


(1 - v)2(l - v(l - w)) 

[(l - 2v(l - w) H- v2(l -2w + 2w2)^ X 
(2 - 2v(l -w) + y2(l - w)2) j, 

1 




- 4w- 


(1 - V)2 


(Al) 




1 


-[2 — w — 2y(l + 4w) 

I L 


^2 


^q^q 

U 3 


^q^q 

'-4 


^q^q 

*^5 


(1 - y)2(l - yvy)(l - y(l - w)) ' 

-i-y2(2 -H 9vy - lOw^ -h w^) - 2y2(l -w + w^- 4w^) 
+v^w{2 -2w - Iw^ + 8 w 2 ) 

-2y^w2(l - 3w H- 4vy2 - 2vy2)j, 

2(1 H-y2(l -h2w2 )) 

(l-v )2 ’ 

-2vw(3 - 2y(l - w) -H y2(l - 2w H- 2vy2)) 

(1 - V )2 ’ 

-T- [2 — w — 2 y( 1 + 2w) 

(1 - y)2(l - vw)(l - y(l - w)) L 

-i-y2(2 -H 5vy - 6w^ + w^) - 2v^(l - w + w^ — 2w^) 

+v^w{2 — 2w- 3vy2 + 4vy2) 

-2yV^(l - 3w H- 4vy^ - 2vy^)], 

1 


, - [2 — w — 2y(2 — w) 

(1 - y)2(l - yw)(l - y(l - w)) ^ 

-i-y2(4 - w - 2w2 + vy2) - 2v^ -b v‘*w(2 - w^) 

-2y^w2(l - 2w-b 2vy2 - w2)j. (A2) 
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'-'1 


'-'2 


r”? 

'-'3 


b. q ^ g channel: 

2vw(l + vHl - wf) 

(1 - v)2(l - v(l - w))2 
X (1 - 2v(l - w) + v^(l - 2w(l - w))), 

vw(6 - 4vw + 2v^(l - 2w(l - w))) 


(1 - V )2 

vw 


r[3-2v(3 + w) 


(1 - v)^(l - vw)^(l - v(l - w))^ * 

+v^(6 + 4w - w^) - 2v^(3 - 3w + 5w^ - 2w^) 
+v‘*(3 - 4w + 5w^ - 2w^) 

-2v^w(2 - 6w + 9w^ - 7w^ + 21 ^“*) 

+2v®(l - w)^w^(l - 2w + 2w^)j, 


r”? 

'-'4 


2 - 2v(5 - 3w) 


(1 - v)2(l - vw)2(l - v(l - w))2 
+v^(16 - 3w - 1 Iw^) - v^(10 + 15w - 27w^ + 2w^) 


+v\2 + 17w - 23w^ 


■ 7w^ 


3w^) 


*^1 


_ 


-v^w(5 - 5w - + 3w^ - 2w‘^) 

+2v®(l - w)^w^(l - w + w^)j. 

g ^ q channel: 

2(1 + v(4vw^ - 2w(l + v) + v)) 


(A3) 


(1 - v )2 


'-'2 


1 


(1 - v)2(l - VW)2 


[2(1 - W + W^) 


'-'3 


-2vw(3 - 2w + 2 w2) 

+y2(2 - 4w + 1 Iw^ - 2w^ + 2w‘^) 

-4v^w(l - 2w + 3w^) + 3v^w^(l -2w + 2w^)j, 

-2v(l +3w + w^ - 6w^) + 


v^(l + 9w + 4w^ ■ 


■8w" 


■ 6w‘^) 


-v^w(3 +9w - 4w^ - 6w^) 
+v‘*w^(l + 4w - 4w^)j. 


(A4) 
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